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Abstract 

We study the differential conductance, spectral density and magnetization, 
for a quantum dot coupled to two conducting leads as a function of bias 
voltage Vds, magnetic field B and temperature T. The system is modeled with 
the Anderson model solved using a spin dependent interpolative perturbative 
approximation in the Coulomb repulsion U that conserves the current. For 
large enough magnetic field, the differential conductance as a function of 
bias voltage shows split peaks. This splitting is larger than the corresponding 
splitting in the spectral density of states pd(w), m agreement with experiment. 
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I. INTRODUCTION 



Electron transport through nanoscale quantum dots (QD's) has been a subject of great 
interest in the last years. A QD consists of a confined droplet of electrons and it has 
been predicted that the system behaves as a single-electron transistor described by the 
Anderson model. [1,2] Experiments with one QD in the linear response regime [3-6] displayed 
clearly the Kondo physics in this almost ideal "one-impurity" system and confirmed the 
predictions based on the Anderson model. [1,2,7] The unitary limit of ideal transmittance 
has been reached. [6] Calculations using the accurate Wilson's numerical renormalization 
group (NRG) have shown a good agreement with these experiments. [8] The effects of a 
magnetic field and the possibility to use transport through a QD or a site coupled to a QD 
as a spin filter have also been addressed theoretically, [9-12] including NRG results. [10] 

In contrast to the linear response regime, the experimental situation in which a finite 
bias voltage between drain and source Vd s is applied, [13-18] i.e. the nonequilibrium case, 
is much harder to treat theoretically. The most accurate techniques used to study the im- 
purity Anderson and Kondo models at equilibrium, [19] cannot be easily extended to the 
non-equilibrium situation. The NRG has not been implemented out of equilibrium. A 
formalism to generalize exact Bethe ansatz results for the nonequilibrium case has been 
proposed recently, [20] but its application to the Kondo model requires further elaboration. 
Some results based on integrabilty but not exact have been presented. [21] Also a Kondo 
model has been solved exactly for particular parameters. [22] The non-crossing approxi- 
mation (NCA), which has been successfully used to calculate thermodynamic properties of 
generalized Anderson models above or near the Kondo temperature T K , [23-25] has also 
been used for nonequilibrium situations. [26-28] Unfortunately, at low temperatures this 
method fails to fulfill Fermi liquid relations in the equilibrium case, [29] and the spectral 
density presents a spurious peak at the Fermi level in the presence of a magnetic field B. 
[27] A likely cause of this peak is given by Moore and Wen, who study the splitting of the 
Kondo resonance with B from the Bethe-ansatz solution at equilibrium. [30] Slave-boson 
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mean-field approximations have been used successfully in the equilibrium case, [12,31-35] 
but they rely on the minimization of the free energy, which is not at a minimum in the 
nonequilibrium case. Nevertheless, these approximations have been used for problems out 
of equilibrium. [36,37] A relatively simple but accurate approach based on perturbation the- 
ory is the local moment approach, [38-40] but its extension to the nonequilibrium case is 
not trivial, because it uses Fermi liquid relations valid at equilibrium. For the Kondo model, 
a perturbative approach in the exchange constant has been used, which is valid when either 
Vds or B are much larger than T K . [41] 

Another approximation for the Anderson model is perturbation theory in U / (vrA), where 
U is the Coulomb repulsion and A is the resonant level width. [42-45] It can be extended 
naturally to the nonequilibrium case using the Keldysh formalism. [46-49] In the equilib- 
rium case, the second order approximation has been used for several nanoscopic systems. 
[12,35,50-53] A shortcoming of this method is that it cannot describe the exponential de- 
crease of T K in the Kondo regime as U is increased. One way to avoid this problem is to 
use renormalized perturbation theory (RPT). [54,55] Another method to extend the validity 
of the approximation to larger values of U is to use an interpolative perturbative approach 
(IPA) , [56-59] which corrects the second-order result in order to reproduce exactly the atomic 
limit U/A — > +oo. The IPA has been shown to describe well the conductance through a 
QD for U < 8 A. [51] The results agree with those obtained more recently using the finite 
temperature density matrix renormalization group method. [60] In addition, comparison of 
the spin dependent IPA [12,52] with exact diagonalization in finite systems shows very good 
agreement for U = 6.25A. [52] Note that in some experimental situations A ~ 0.15 meV, 
U ~ 0.6 meV (U/A = 4), [3,5] while in recent nonequilibrium experiments 2 A ~ 0.33 meV, 
U ~ 1.2 meV (U/A = 7.3). [16] The extension of perturbation theory in U 2 to the nonequi- 
librium case has been considered by Hershfield, Davies and Wilkins. [61] They found that 
for finite bias voltage Vd s , the current is conserved only in the symmetric Anderson model, 
in which the on-site energy at the dot E d = —U/2 + (/i L + Hr)/2, where /il (/Ar) is the 
chemical potential at the left (right) lead. Recent use of perturbation theory in U / (nA) for 
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Vds 7^ was limited to the symmetric case, for calculations of the current noise, [62] and 
other properties including terms of third and fourth order in U, [63-66] and in particular 
in presence of magnetic field B, [64,65] motivated by recent experiments. [14-16] However, 
even in the symmetric case, the spin current is not conserved by the approximation when 
Vd s 7^ 0. In the Kondo regime, an exact expression for the differential conductance has been 
derived using RPT, [65] but it is limited to Vd s small in comparison with Tk- 

There is some controversy concerning the effect of Vd s on the Kondo peak in the spectral 
density for B = 0. Calculations based on the NCA [26-28] and on the equation of motion 
method [67,68] predict a splitting of the peak. However, this method has limitations in 
the Kondo regime, particularly in the particle-hole symmetric case, for which the Green's 
function is indepedent of temperature. [69] Instead, second-order perturbation theory in U 
predicts a fading of the peak with increasing Vd s without splitting. It has been suggested that 
terms of order U 3 and U 4 might cause the peak to split. [63,64] However, a recent study with 
an improvement of the higher order corrections, which reproduces correctly the equilibrium 
expressions, indicates that the spectral density retains the qualitative features of the second- 
order result. [66] Physically, one might expect that if the dot is hybridized with left and right 
leads at chemical potentials /il and jiR with matrix elements Vl and Vr respectively, the 
effective distribution function fluctuates between that corresponding to hl and /ir (at least if 
Vl, V r > \/jLl — Hr\) and no split Kondo peak at the Fermi level is expected. [66] However, as 
stated above, different approximations lead to different results. The experimental situation is 
also controversial: a fading of the central peak in the differential conductance with increasing 
bias voltage was reported for a carbon nanotube quantum dot, [17] while a splitting of the 
Kondo resonance in the spectral density was observed in a three-terminal quantum ring. [18] 
Since the diameter of the ring is of the order of 500 nm, it might be possible that to model 
the system as one effective site connected to conducting leads is not a good approximation 
and the space dependence of the energy distribution function plays an essential role. [70] 

In this work, we generalize the spin dependent IPA, based on perturbation theory in 
U up to second order, to the nonequilibrium case. The problem of the conservation of 
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the current for each spin is solved by a simple trick. The partition of the Hamiltonian 
H — Hq + H' into an unperturbed part H and a perturbation H' is in principle arbitrary. 
In the past, for example H has been chosen so that the Friedel-Langreth sum rule [71] 
is satisfied at zero temperature. [51,59] Here we also use the spin dependent version of 
this rule [12,54] to calculate the spectral density at equilibrium. For finite bias voltage 
Vds, we determine H asking that the current is conserved for each spin projection. The 
formalism is applied to calculate the differential conductance dI/dVd s , spectral density and 
magnetization as a function of bias voltage and magnetic field B for various temperatures. 
Our results are compared with experiment and previous calculations. Recent measurements 
[14-16] of dI/dV ds report results which seem to disagree with accurate calculations [72,73] 
of the spectral density at equilibrium for B ^ 0. However, Hewson et al. pointed out that 
these experiments should not be interpreted in terms of the equilibrium Green's functions. 
[65] In spite of shortcomings of the selfconsistent procedure at small temperatures, when 
both Vds and B are small but different from zero, important conclusions can be drawn from 
our results. 

The paper is organized as follows. In Section II the model is described, the nonequi- 
librium perturbation formalism is briefly reviewed, useful expressions are derived, and the 
generalization of the IPA to the nonequilibrium situation is presented. Section III contains 
the results of the application of the method to magnetotransport, spectral density and mag- 
netization. Section IV contains a summary and a brief discussion. Some details of the 
calculations were moved to the appendix. 

II. MODEL AND APPROXIMATIONS 
A. Hamiltonian 

We consider a QD interacting with two conducting leads, one at the left and one at the 
right, with chemical potentials hl and /xr respectively, with hl-^r = eVds- As usual, the QD 
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is modelled by one effective site with one relevant orbital with an important on-site Coulomb 
repulsion U and an on-site energy E d controlled by the gate voltage V g . The Hamiltonian is 
that of an Anderson model, which is split into a noninteracting part H and a perturbation 
H' as: 

H = H + H', 

kaa a kaa 

H' = (E d - £eff)^2 n da ~ B(n d ^ - n dl ) + Un d ^n dl . (1) 

a 

Here a = L,R refers to the left and right leads. The operator c\ aa creates an electron in 
the state with wave vector k and spin a at the lead a. Similarly d\ creates an electron 
with spin a at the QD. The number operator n da = S a d a , and B is the effect of an applied 
magnetic field on the on-site energy at the QD for each spin The first term of H describes 
the leads and the third its hybridization with the QD. H' contains the dot on-site energy, 
Zeeman splitting and Coulomb repulsion respectively. The two constants e c e ^ are arbitrary 
since they cancel in H + H' . However, in our perturbation treatment, the results depend 
on them and they are determined selfconsistently as described in subsection E. 

B. Nonequilibrium perturbation theory 

We use the notation of Lifshitz and Pitaevskii. [47] There are four different Green's 
functions. For the electrons at the dot, they can be chosen as 

GL(M') = ~i9(t - t')(d a (t)dUt') + 4(t')da(t)), 
G a da (t,t') = iO(t - t)(d a (t)dt(t') + 4(<%(*)>, 
Gl(t,t')= i (dUt')da(t)), 

G> a (t,t>) = -i(d a (t)dt(t')). (2) 

The first two are the retarded and advanced Green's functions, already present in the equi- 
librium case. Similar Green's functions can be constructed involving conduction electrons. 
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Since obviously G a — G r = G K — G > , only three Green's functions are independent, and 
in the following we use this relation to eliminate G > . Also, three independent self-energy 
functions S r , S a , and can de defined, which allow us to write a matrix Dyson's equation 



[47] 



G = g + gSG, 



(3) 



where 



G = 



/ 



V 



G a 
G r 2G < + G 1 



Qa 



(4) 



g has a similar expression in terms of the corresponding Green's functions for the unper- 
turbed Hamiltonian H , and 



V 







(5) 



/ 



In general, the products in Eq. (3) have to be understood as convolutions in time and 
space. However, in this work, in which we consider the stationary case and perturbation 
theory in U up to second order, [66] these products become just ordinary products involving 
the Fourier transforms in t — t' of the Green's functions for the d electrons G v da {uj) and the 
corresponding self energies In addition, G a da {uj) = (G r da (uj))* and S"(cj) = (S^(cj))*, 

where the asterisk denotes complex conjugation. Therefore, our task reduces to find suitable 
approximations for S^.(cj) and S^(cj). 

After doing the matrix product in Eq. (3), from the entries (1,2) and (2,1) one obtains 



-i 



(6) 



and its complex conjugate. After some algebra and using previous equations, the remaining 
equation (from the entry (2,2)) can be written in the form 



g < m . |g ^ )|2 (^_ e;m ). 
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(7) 



The noninteracting Green's functions can be obtained easily from the equations of motion 
[49] using H : 

a 

where f(u) is the Fermi function and E a , A a for leads a = L,R, are the real and imaginary 
parts of the sum 

V |Vte| * . = E a {u) - iA a (u;), (9) 
u - e ka + ir] 

with i] a positive infinitesimal. 

We approximate the retarded self energy as 

£» = El - e° eff + U(n^) + (10) 

where E\ = E d - B, E l d = E d + 5, and a means —a. The last term is an interpolative 
expression based on the correction of order U 2 , as described in subsection D. The remaining 
terms would correspond to the first order correction in H' if (n^) were evaluated with H . 
However, (n^) is evaluated selfconsistently using 

M = duG<M, (ii) 

and this is equivalent to a partial summation of an infinite series of diagrams. [64] 

The correction for £<(cj) of first order in H' vanishes. We approximate S<(c<j) by an 

expression £< int (u;) based on the second order term E^cu), also described in subsection D. 
The diagrams for the corrections to the self energy of second order in U are drawn in 

Ref. [47]. In terms of the independent unperturbed Green's functions the expressions are 

[g^MglMgfi^i + u 2 - u) 
+g\MgfMgf{^i + u 2 -uj) 
+gfMg r i Mgf(ui + u 2 -uj) 

+g<(u l )g<(uj 2 )g a l (u l +u 2 - u)], (12) 



E T <2 (cu) = -U 2 



du>i f du- 



2tt J 2tt 

gf(u)i)gf(u 2 )gl(uji +u 2 -w), (13) 

and the same interchanging spin up and down, with g^{ui) = g^{oj) + 2lmg^(uj). It can be 
shown that these expressions are equivalent to those given by Hershfield et al. [61] These 
integrals which have the form of convolutions in frequency can be calculated conveniently 
using fast Fourier transform to the time representation. [45,74] For the case of constant 
E a (u) and A a (u), it is more convenient to evaluate one of the integrals analytically, as 
shown in the appendix. 

Although it is not convenient for the numerical evaluation, it is interesting to express 
the above self-energy corrections in terms of the unperturbed spectral density of states at 
the dot Pda(ou) = — \va g r du (uj) / it . In terms of this density 



9 r M = [ de-^ 

J LU — t 



e + ir\ 

5 <H=2i7rpLH/», (14) 
where f(cu) is a weighted distribution function of the Fermi functions at the two leads: 

/» = s« A °<^-^ (15) 

Replacing Eqs. (14) in Eqs. (12) and (13) we obtain after some algebra 

J uj + e 3 - ei - e 2 + %r\ 

x[(l-/(6l))(l-/(6 2 ))/(63) 

+/(ei)/(e 2 )(l-/(e 3 ))], (16) 

£<» = -2tnU 2 J de 1 de 2 p%(e 1 )p dl (e 2 )p° dl (e 1 + e 2 - u) 

x/(ei)/(e 2 )(l-/(ei + e 2 -a;)). (17) 

At equilibrium pl — Pr- Then f{uS) becomes the usual distribution function and Y7 2 (uj) 
coincides with the result of ordinary perturbation theory. Eqs. (16) and (17) are used in 
subsection D to construct the corresponding interpolating expressions E£f nt (u;) and £< int (cj). 
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C. The current 

Following Meir and Wingreen, [75] the current with spin a flowing between the left lead 
and the dot is 

lie f 

ha = T J duA L (u) [2if(u - p L ) ImG^» + G%(u)\ . (18) 

Similarly, the current with spin a flowing between the dot and the right lead is 

lie f 

JR. = — J cLuA R (cu) [2if(u - i* R ) Im(7» + G% (u)] . (19) 
Of course, since the current is conserved one should have 

jL*=jRa, (20) 

but perturbation theory does not satisfy this equality for generic values of the parameters. 
[61] In this work we determine e a e ^ imposing Eqs. (20). 

D. The interpolative perturbative approach (IPA) 

Extending previous ideas, [56-59] we want to replace the second-order contributions to 
the self energies Y7^{uj) and £< 2 (a>), by other ones ££.f nt (a;) and E^ int (o;), which coincide 
with the previous ones to order U 2 for small U, but also reproduce the high frequency and 
atomic limits. For the sake of clarity in the exposition we choose a =f. The final results 
are also valid interchanging spin up and down. 

We propose [59] 

E «<-> - (21) 

where A 1 is determined so that £|f nt (u;) reproduces the leading behavior at high frequen- 
cies, and afterwards A 2 is determined to reproduce the exact result in the atomic limit 
A a (u)/U^0. 
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Up to order 1/cu, the exact self energy is determined by the first and second moments of 
the spectral density of d states p d(T (u), which can be evaluated independently of p da {uj) using 
particular commutators. [76] Proceeding as in the equilibrium case, [59] allowing dependence 
on spin, [12,52] and using Eqs. (6) and (10), one obtains that the correct leading behavior 
of the retarded second-order self energy is 

Eg» = (n dl )(l - (n*))!/ 2 *" 1 + 0(c- 2 ). (22) 

For uj — > oo, the integrals over e, in Eq. (16) decouple and then 

Sf H = (n%)(l - {n di ))U^ + 0(uj- 2 ), (23) 

where (n da ) = (—i/2n) f dcug^cu) is the expectation value of n da calculated using the 
noninteracting lesser Green's function [see Eq. (14)]. From the above equations, one sees 
that choosing 

A 1 = (n dl )(l-(n d J)/[(n dl )(l-(n dl )}, (24) 

the moments of p du {u) up to the second one are reproduced exactly. 

In the atomic limit, it can be easily verified using equations of motion, [49] that the exact 
retarded Green's function is 

G d f{u) = 1 = ^ + . (25) 

T cu-El + ir) cu-Ej-U + ir) 

From Eqs. (6) and (10), this implies 



u;-El-(l-{n dl ))U + irj 



^Tint K U ) ~ pf 7, ; u„ . . • \ Zb > 



The usual contribution in U 2 to the retarded self energy can be calculated easily in the 
atomic limit using Eq. (16), because the noninteracting spectral densities p da {oj) become 
delta functions. The result can be written in the form 

Ef -V) = (27) 
uj - e[ ff + ir] 
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From Eqs. (21), (24), (26) and (27), the coefficient A 2 is determined. The final result is 



T (<)) + [4- £? l-( 1 -WMW/y J )' 

This result has the same form as in the spin dependent situation in equilibrium. [12,52] 
However, the expectation values are calculated with distribution functions out of equilibrium. 

From the equations of motion in the atomic limit (see section 3 of Ref. [49]), one obtains 
that the lesser Green's function for Vk a = is the sum of two delta functions G^ at (uj) = 
iaS(uj — E\) + ib5(u — E\ — U), where a and b are two real unknown constants. Using Eq. 
(25), we can write G^' at (cj) = — 2iy(uj) IrnGj^w), where the function y(u) is unknown. 
Electron-hole symmetry in the symmetric Anderson model imposes that y{—uo) = 1 — y(ou). 
At equilibrium, \i L = /i R = \i and it can be shown from their definitions [Eqs. (2)] [75] 
that y{oS) = f{oJ — /i). Out of equilibrium, y{uj) is ill defined at V^ a = 0, since the dot is 
disconnected to the leads, while even for small non-zero Vk a the problem is not trivial. In 
the general case, we assume that y{oj) = f(u>), the average distribution function defined in 
Eq.(15). This has the correct limit at equilibrium or when the limit Vk a — > is taken for 
one of the leads first and then for the other. Using this choice, Eqs. (6), (7) and (8) imply 
for the lesser self energy in the atomic limit 

S< n tH = 2./»ImSgfH. (29) 

On the other hand, using Eqs. (8), (17) and (27), we can write for the second order correction 
in the same limit 

S^' at (u) = 2if(u) Im Sf ' at (u) . (30) 

From these equations, and taking into account that E|? nt (u;) coincides with Y^{u) for small 
U and with S^ 2 ' at (cj), the simplest interpolative expression for S^(cj) is 

ImEj.(w) 

^TintM = ImSfM S T (31) 

when Im S| 2 (w) 7^ 0, and zero otherwise. Note that if ImS| 2 (w) = 0, also E| < (c<j) = 0, as 
can be seen from Eqs. (16) and (17). 
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E. The selfconsistency problem 



Eq. (28) replaced in Eq. (10) and the same interchanging j and J, define the retarded self 
energies. Similarly Eq. (31), and the same for spin down are the lesser self energies used. 
The expressions for Yi^iuj) and T^{uj) for constant A Q are given in the appendix. The self 
energies replaced in Eqs. (6) to (8) give the Green's functions for the interacting system 
that depend in general on four unknown quantities: e a e ^ and (n da ). They are determined 
by the selfconsistent solution of Eqs. (11) and (20). The resulting values of e a e ^ do not have 
any particular physical meaning. 

In the present work, we present results for the symmetric Anderson model, E d = -U/2, 
ix l = — (J>r = eV ds /2, Al = Ar independent of energy, in presence of a magnetic field B. 
In this case, due to electron-hole symmetry e a e ^ = — ££// an d 3 La — 3m- This reduces the 
problem to two selfconsistent equations with two unknowns. 

Unfortunately, for small temperature T, small non-zero V ds and small non-zero B, there 
were instabilitities in the numerical algorithm used. Also in some regions no solution could 
be found, while in other cases the solution jumped to another value. Therefore our results 
should be regarded as complementary to those obtained using renormalized perturbation 
theory (RPT), which correspond to T = and small V ds . [65] We have tried an alternative 
approach, relaxing the condition of conservation of current and fixing e a e ^ either at the value 
that satisfies the Friedel-Langreth sum rule (see below) or e c e ^ = 0. In the former case, the 
solution of Eqs. (11) for (n da ) was lost for very small V ds , while in the second, the peak 
in the differential conductance dI/dV ds was already split for very small values of B (of the 
order of TV/30) in contradiction with the results of RPT. Therefore these alternatives were 
abandoned. 

In general, we have started from the solution for large V ds (with small reasonable ke//|), 
which was easy to find, and used it as an initial guess for lower V ds . The procedure was 
repeated until reaching a small non-zero value of V ds or some instability. 

For V ds = 0, the problem is ill posed since the current is conserved for any ef^. Also 
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the first derivative of jha — 3ro with respect to V ds vanishes and the second is very small for 
any e a e ^. For V ds = T = 0, we have calculated the spectral density replacing the condition 
jia = ]R a by the Friedel-Langreth sum rule [71] for non-zero magnetic field, which for 
constant A = A R (oo) + A l (lu) can be written in the form [12,54] 

nAp df7 (0) = A = sin 2 (7r(n, CT )). (32) 

A2 + Kf + ^M] 



III. NUMERICAL RESULTS 

In this section we take A R (u) = A L (u) independent of energy and A = A^ + A L = 1 
as the unit of energy The origin of energy (cu = 0) is set at (pl + /j._r)/2. We have chosen 
U = 5 for the numerical analysis because we want U large enough so that the system is at the 
Kondo regime, but for U > 6.5, the magnetic susceptibility for V ds = T = decreases with 
increasing U, indicating the failure of the IPA for a quantitative analysis at small non-zero 
magnetic field. These parameters lead to a Kondo temperature T K = 0.45 as determined 
from the width of the spectral density. 

To represent the spectral density for V ds = T = 0, we have determined selfconsistently the 
occupation numbers (n da ) and effective energies e u e ^ using the Friedel-Langreth sum rule, 
while to calculate the conductivity at finite V ds we have replaced this rule by the condition 
of conservation of current as described above. 

In Fig. 1 we represent the differential conductance dI/dV ds as a function of the applied 
voltage between drain and source V ds for zero temperature, and compare it with the total 
spectral density of states p d {uj) = p d -\(uj) + p d i(uj) for V ds = and uj = eV ds . The spectral 
density at B — shows a central peak corresponding to the Kondo peak with a half width at 
half maximum T K = 0.45 and two shoulders at uo = E d = —2.5 and u — E d + U = 2.5, which 
tend to separate with increasing magnetic field. The central peak shows a splitting which 
is tiny for B = 0.2 (not shown) and increases with applied magnetic field. The splitting of 
the peak is larger than the Zeeman splitting 2B by a factor ~ 1.5 for B — 0.3 and increases 
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to slightly below 2 (~ 1.8) for B — 1.5. This behavior agrees with previous studies of the 
spectral density under an applied magnetic field for large enough U [10,26,30,39,55,64] 

The differential conductance dI/dVd s was obtained by numerical differentiation of the 
total current / = 2l\ + 3l{ = jni + Jri- For < B < 0.5 and < V ds < 1-5, we could not 
find a self consistent solution of the system of equations, and therefore we are not able to 
describe how the peak in the differential conductance splits with applied magnetic field at 
T = 0. As described below, the situation improves with applied temperature. For B = 0, 
the peak in dI/dVd s is narrower than that of the spectral density with a half width at half 
maximum of 0.383. In addition the shoulders are displaced from the central peak. In general, 
the split peaks in dI/dV ds are displaced to higher \eVd s \ when compared with the same peaks 
in pdiyj). Because of the lack of results for small non vanishing magnetic field and voltage, 
our calculations for T = are not enough to establish if this tendency persists for small 
splittings. 

In Fig. 2, we show the evolution of the spectral density with applied voltage. For B = 0, 
the Kondo peak decreases slowly without splitting. For Vd s = 3 (not shown), Pd{oj) is flat 
near u = 0. Further increase in Vd s leads to only very small changes in p d (u). For large 
magnetic field, the changes in Pd{^) with V ds are less dramatic. For large V ds the spectral 
density for different applied magnetic fields looks similar, but somewhat flatter and broader 
for larger B. 

In Fig. 3 we show the magnetization m = (n^) — (n d i) as a function of applied voltage 
between drain and source for different values of the magnetic field and temperature. Because 
of our problems with the self-consistent solution for T = 0, the values of B > T K = 0.45 taken 
at zero temperature are large enough so that the Kondo effect is at least partially destroyed. 
In this case, an increase in Vd s leads to a decrease in the magnetization which is expected as a 
consequence of the flattening of the spectral density of states discussed above. These curves 
for T = and comparatively large values of B show an initial slow decrease with increasing 
applied voltage, an inflection point for eVd s ~ 45, and an asymptotic approach to m = 
for large Vd s - Instead for small B (see the curve for B = 0.2 and T = 3), the magnetization 
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first increases and then decreases with applied voltage. This can be interpreted in terms of 
a destruction of the Kondo singlet state for V ds < 2T K , which moves the system towards 
the local moment regime, followed by the same effects of flattening of the spectral density 
present for higher B. 

For temperatures T ~ T K /2, we have obtained some results in the problematic region 
for selfconsistency (0 < B < Tk and < Vd s < 27V). Some of them are shown in Fig. 4. 
However, still the value of B at which the splitting of the central peak takes place could 
not be identified precisely and abrupt changes in the solution take place for example for 
B = T = 0.3 and B = T = 0.2. In spite of this, these results and those for T = 0.5 shown 
below, suggest a transition from one maxima to two maxima in the differential conductance 
dI/dV ds as a function of applied voltage for B c ~ 0.3 = (2/3)T K . This value is consistent 
with that obtained by Hewson and coworkers [65] B c = 0.5847^. Comparison of all split 
maxima of dI/dVd s with the corresponding ones of the spectral density, shows again that 
all the former lie at values of \eVd s \ larger than the corresponding values of \u\ for the 
maxima of Pd(yj). This provides a natural explanation of the disagreement of the splitting of 
the peaks observed in experiments of the differential conductance in presence of a magnetic 
field, [14-16] when compared with accurate results of the spectral density [72]. The observed 
splitting for sufficiently large field was found to be higher than the corresponding one of the 
spectral density. 

In Fig. 5, we show the evolution of Pd{u) with applied voltage for T = 0.3. For B = 0.2 
and small V^ s , the spectral density displays only one maximum at w = 0, although its 
components p d -\(uj) and p d [(iv) have maxima at uj = —ujm and uj = % respectively, with 
ujm > 0, as shown in Fig. 5 (a) for eV ds = 0.1. We remind the reader that due to our choice 
of parameters Pdi(^) = Pd](— w). As in the case of zero temperature (Fig. 2), the changes 
driven by the applied voltage on Pd a {u) are smaller for higher magnetic fields. 

In Fig. 6 we show dI/dV ds as a function of Vd s for T = 0.5, and the spectral density pd{oS) 
for eV ds = 0.1. As before and to facilitate comparison, the abscissa of p d (uj) is chosen in such 
a way that uj in Fig. 6 (a) is the same as eVd s in Fig. 6 (b). The scales are the same as in 
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Fig. 4. For this temperature, there were no instabilities in the algorithm for selfconsistent 
equations and we could follow the splitting of the central Kondo peak. The transition from 
one maxima to two in the differential conductance takes place slightly below B = 0.3. With 
increasing temperature, both the spectral density and the differential conductance become 
flatter, and the changes are more noticeable for small vales of the magnetic field. Again, the 
structure of dI/dVd s is broader than that of Pd(a>) and the maxima are more separated from 

LO = 0. 

IV. SUMMARY AND DISCUSSION 

We have generalized the spin dependent interpolative perturbative approximation (IPA), 
based on second-order perturbation theory in the Coulomb repulsion U, to the nonequilib- 
rium case. This permits to extend the validity of the results of nonequilibrium second order 
perturbation theory in U / (VA) to higher values of U. The effective unperturbed energy at 
the dot has been determined selfconsistently asking that the current for each spin is con- 
served. This allows us to correct the shortcoming of ordinary perturbation theory in U/ (7rA) 
that the spin current is not conserved, except for the symmetric Anderson model without 
an applied magnetic field. 

We have applied the method to study the effects of a magnetic field B on the differential 
conductance dI/dVd. s , starting from the symmetric Anderson model in which the on-site 
energy E d = —U/2, setting the origin of energies at the average of the chemical potentials, 
and the resonant level width for both leads are equal (Al = Ar) and independent of energy. 
The selfconsistent procedure fails for temperatures below the Kondo temperature T K , when 
< B < T K and < V ds < 2T K . 

The procedure can in principle be applied to the general case (without starting from 
the symmetric case). Then, one has to solve four self consistent equations for the effective 
occupations up and down, and the corrresponding effective unperturbed levels. Similar 
problems were solved for the equilibrium case in presence of a magnetic field [12] and to 
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calculate persistent currennts in a ring with an odd number N of electrons. [52] In the latter 
case however, the procedure failed for large N. In the present case, one might expect more 
numerical difficulties after our experience with the simpler case treated here. 

In spite of the above mentioned failure, our results indicate that the critical magnetic 
field B c for the splitting of the Kondo resonance in dI/dV ds as a function of V ds is near 
and slightly below (2/3)7^. This is consistent with the results of renormalized perturbation 
theory [65] which give B c = 0.5847V. 

In the region of split peaks, the splitting is larger than the corresponding one of the 
spectral density of states paiyj) as a function of uj. This provides a natural explanation 
of why recent experiments of the differential conductance in presence of a large magnetic 
field give a splitting larger than two times the Zeeman splitting, as expected for the spectral 
density of states of states at equilibrium. [14-16] Instead, previous results using slave bosons 
obtained a reduction of the peak splitting. [37] 

The effect of increasing bias voltage on the spectral density Pd(w) is just an overall 
flattening of the energy dependence, without introducing new peaks. This is in agreement 
with recent calculations in the symmetric case [66] including terms of third and fourth order 
in U. 

For small applied magnetic field, so that the system is in the Kondo regime, the mag- 
netization as a function of the applied bias voltage V ds first increases because due to the 
destruction of the Kondo singlet state, the system enters in the localized moment regime. 
For Vd s ^ 27V, or B > TV, the effect of flattening of the spectral density dominates and the 
magnetization decreases with increasing Vd s - 
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APPENDIX A: SELF ENERGIES FOR CONSTANT A, 



Usually for problems in QD's, the density of states of the wide band of conduction 
electrons at the leads, as well as the hybridization with the dot Vk a can be taken as constant 
within the narrow energy range of Tk or other energy scales in the problem. [51] In this case 
E a and A a become independent of energy and one of the integrals in Eqs. (12) and (13) 
can be done analytically E a merely renormalizes Ed- 

Decomposing g^ioj) as well as products of unperturbed Green functions as sums of 
simple fractions with single poles times Fermi functions, and decomposing products of Fermi 
functions using 

! ^= ^ [(X+ \),T]-^ -^ < A1) 
the first integral in frequency is decomposed into a sum of terms of the form 

duf(u) (A2) 

uj — a ± ib 

with a and b real. This integral can be evaluated using the digamma function: [77] 

/ °-/^)^=*4 + w)-*'(^ (A3 » 

After a lengthy algebra, the corrections of order U 2 to the self energies become 

+ g <(u;-v)h(v)l (A4) 



2fM = -(^) / dug<(u-u)h(u), 



and the same interchanging spins up and down, with 

2 „ , 2iA 



h(u) = — Re{ ^ + 2?A) J2U a + B & W 2a {-u) - if**) (A5) 



+(A a + B a )(<p 2a (v)-<p* la )}}, 

2iA a 
v(v + 2iAy 
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where A = A^ + A^, a = R (L) if a = L (R), and ipi a and the functions tp2a(v), A a {y) and 
S a (z/) are 

1 A + ifc^-^J 
^ = ^2 + 2^T } ' 

,1 M^-^ 

^W=^2 + £t } ' 

A 2 

A» = 



exp(z//T) - 1' 



= F7 Al/AR r ; ; . (A6) 

In equilibrium, when Hl — I^r, the expression for YT^(uo) can be easily shown to coincide 
with that given by Horvatic and Zlatic. [43] 
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FIGURES 




FIG. 1. (a) Spectral density at the quantum dot as a function of energy for Vd s = and (b) 
Differential conductance as a function of applied voltage, for T = 0, and different values of B. We 
remind the reader that Tk = 0.45. 
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FIG. 2. Spectral density as a function of energy for T = 0, different values of the applied 
voltage V = eV ds and (a) B = 0, (b) B = 1. Note that Tr- = 0.45. 
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FIG. 3. Magnetization as a function of the applied voltage for T = (except dashed dot line 
for which T = 2T K /3 = 0.3) and several values of B. 
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FIG. 4. Spectral density a function of energy for Vd s = 0.1 and (b) Differential conductance as 
a function of applied voltage for T = 27V/3 = 0.3, and different values of B. 
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FIG. 5. Spectral density as a function of energy for T = 27V/3 = 0.3, different values of 
the applied voltage V = eV^s and (a) B = 0.2, (b) B = 0.4. The thin line corresponds to the 
contribution with spin up pd\(u) for V = 0.1. 
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FIG. 6. Same as Fig. 4 for T = 0.5 ~ T K = 0.45. 



31 



